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Abstract 

We assess the model discriminating power of a combined phenomenological anal- 
ysis of /X — > 67 and /i — > e conversion on different target nuclei, including the current 
hadronic uncertainties. We find that the theoretical uncertainties can be largely re- 
duced by using input from lattice QCD and do not constitute a limiting factor in 
discriminating models where one or at most two underlying operators (dipole, scalar, 
vector) provide the dominant source of lepton flavor violation. Our results show that 
a realistic discrimination among underlying mechanisms requires a measurement of 
the ratio of conversion rates at the 5% level (two light nuclei) or at the 20% level 
(one light and one heavy nucleus). We have illustrated these main conclusions also 
in the context of a supersymmetric model. 



1 Introduction 



Lepton Flavor Violating (LFV) decays of charged leptons provide a theoretically clean 
probe of physics beyond the Standard Model (SM), due to the un-observably small branch- 
ing fractions (~ 10^^°) within the SM, minimally extended to include massive neutrinos. 
Searches for SM forbidden muon processes, such as /i — 67, /i eee, and /i — e con- 
version in nuclei, have provided so far the strongest constraints on LFV dynamics, with 
90% C.L. upper limits given by B^^^e-y < 1-2 x IQ-^^ [I], 5^+_e+e-e+ < 1-0 x lO'^^ [2], 
5^^e(Au) < 8 X 10-^^ [3], 5^^e(Ti) < 4.3 x lO'^^ % 5^^e(Pb) < 4.6 x lO"" [5] 0. 

It is a well known fact that while the decay — > 67 is only sensitive to a transition 
magnetic dipole operator, both fi — > eee and /i — > e conversion in nuclei are sensitive 
to transition charge radii operators as well as purely contact four-fermion interactions 
induced by physics beyond the SM. In other words, different LFV decays have different 
sensitivities to underlying LFV mechanisms (effective operators). This leads naturally to 
ask the question whether one could infer the relative strength of these different operators 
in a completely phenomenological and model-independent way. This would allow one to 
discriminate among different underlying models of LFV and thus would provide valuable 
input for model building. 

In Ref. |S] it was pointed out that in principle, by combining the rates of /i — 67 and 
/X — i> e conversion on different target nuclei, one could obtain information on underlying 
models. There are three types of effective operators that contribute to the coherent fi —>■ e 
conversion process: the dipole, the vector, and the scalar operators. In the non-relativistic 
approximation of the muon wave function, the three operators give the same form of 
overlapping integrals among the wave functions of the initial muon and the final electron 
and the nucleon density in the target nuclei. However, as the relativistic and finite nuclear 
size effects become important for heavy nuclei [H [TJ E] , the transition amplitudes for the 
three operators show different dependences on the atomic number Z. The relative numbers 
of neutrons and protons also change as Z increases. This fact helps to find out if the lepton- 
flavor-violating operators couples to up-type or down-type quarks again by looking at the 
target atom dependence. In this work we go back to this issue with the aim to 

• quantify the theoretical uncertainty induced by the quark scalar density matrix ele- 
ments in the nucleon; 

• quantify the experimental precision required to realistically infer useful information 
on the underlying LFV mechanisms. 

We organize our discussion as follows: in Section[2]we review the derivation of the /i — >■ e 
conversion rate starting from a general effective theory description of the LFV physics. In 
Section [3] we explore the phenomenological consequence of the simplest possible models, 
in which only one effective LFV operator dominates. We extend this analysis in Section 
m to the class of models in which two operators dominate. In Section [5] we specialize 

^ B^^e{Z, A) represents the ratio of — > e conversion rate over muon capture rate, namely ^"""llz'^l ■ 
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our discussion to a supersymmetric (SUSY) model and summarize the conclusions of our 
analysis in Section [61 



LFV effective interaction and the /i 
rate 



e conversion 



In this section we review the procedure to calculate the rate of the fi ^ e conversion in 
nuclei, starting from a general parameterization of new physics effects via effective operators 
at a scale A larger than the electroweak scale v ~ 174 GeV. 



2.1 Effective Lagrangian 

We start with the most general effective Lagrangian which describes LFV transitions be- 
tween charged leptons of first and second families at the weak scale: 



^eff 



1 

"A2 



+ {CGRm^Gp ePifi + CcLm^GF eP^/i) ^G':G% + h.c. 



(1) 



We have not included operators involving q'jp'y^q, ^75^', or qapi,q since they do not con- 
tribute to the coherent conversion processes. In the above expression A represents the scale 
where new physics effects appear. We take A = 1 TeV in this paper. The Cab's are dimen- 
sionless constants containing information about the underlying theory; the subindexes R, L 
correspond to the chirality of the final electron which is determined by Pr^l = (1 =t 7^)/2, 
q are light and heavy quarks. The field strength of the photon and the gluon are defined 
by Fp^ = dpAj, — d^Ap and = dpG^ — d^Gp — fabcG'^pGl, respectively. The normal- 
ization is chosen so that the kinetic terms are given by — (1/4)FF and —{l/4:g1)G°'G°'. 
The a matrix is defined by a^'^ = ^[7^,7^]. Gp = l/(2v^'i^^) is the Fermi constant, while 
rrip and rriq represent the muon and running quark masses at yU = A, respectively. We 
have introduced the running quark masses and the beta function of the QCD coupling 
constant, /? = {g^ / IGn'^) {11 — 2Np/3), so that the coefficients C's do not depend on the 
renormalization scale under QCD running at 1-loop level. The notation (3h,l is used to 
distinguish the Lagrangian with all quarks contributions (H) from the one where heavy 
quarks are integrated out (L). The Lagrangian in Eq. [1] describes three kind of interactions 
that violate the lepton flavor: The effective interaction with a photon (Dipole term), the 
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effective interaction with quarks (Scalar and Vector terms) and the effective interaction 
with gluons (Gluon term). 

In order to evaluate the — >• e conversion rate, it is appropriate to use the effective 
Lagrangian at the nucleon level [9] . We first integrate out the heavy quarks before matching 
to the nucleon level Lagrangian. It can be straightforwardly done by using the matching 
of the trace anomaly [10] . The Lagrangian is given by 

- 
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+ 

q=u,d,s 



=u,d,s 

m^niqGF ePjifi ) qq 

q=u,d,s 

+ {CGQRm.Gp ePLfi + CcQLm.GF ePufi) ^G^^-^G^ + h.c. 



(2) 



where jSi is the beta function of three-flavor QCD. The new coefficients of the gluon terms 
are expressed in terms of the original Lagrangian parameters as follows: 

Cgqr = ^ C^sr '^Q + Cgrh 

Q=c,b,t 

Cgql = ^ 4T '^Q + Cgli^ , (3) 

Q=c,b,t 

where 

_ A{(3/gf) _2_ _ [P/gDn _ 7 

{(3/gl), 27' iP/gl), 9 

The Lagrangian ([2]) can be evolved with the renormalization group down to energy 
scales of the order of /i ~ 1 GeV, by simply taking the quark masses and gauge coupling 
constants in the Lagrangian to be the running ones at ~ 1 GeV. At this low scale, 
we match to the effective Lagrangian written in terms of the relevant degrees of freedom, 
namely nucleons, leptons, and photons. That can be done by the following replacements 
of operators: 

f^q qq fs% i^N^^N 

qipq fvN^NlpipN 

P^GG Jgn ruN ipN'tpN ■ (5) 

where N represents each nucleon (A^ = p,n), ipj^f are the nucleon fields, and /'s are nucleon 
form factors. The form factors depend in principle on the momentum transfer, which we 
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will neglect as it is smaller than the typical scale of the nucleon structure. The fact that 
(A'"|^"|A'") = {Nlt/j^t/jNlN) {92 is the trace of the energy momentum tensor) imphes 
the simple sum-rule 



/ , JSN 
q=u,d,s 



1= V fZ + f. 



GN 



(6) 



which we use to eliminate the form- factor foN in terms of the scalar nucleon form factors 



f^glf. The nucleon vector form factors are known from the vector current conservation, 



Au) 
JVp 



Au) 

JVn 

(d) 



JVp ~ ^ JVn 

f (*) _ n 

JVp — ^ JVn 



1 

2 




(7) 



while the calculation of the scalar form factors f^^ is non-trivial. As discussed below, in 
our analysis we will use input from Chiral Perturbation Theory and the lattice QCD to 
asses the impact of current and future uncertainties on the conversion rate. 

Collecting the above results, the Lagrangian at nucleon level can be written as 

-- y 



^eff 



N=p,n 

W ^.,PP„„ ^^,pu 



ePLi^ + Cf^^ ePnii) V'ivV'iv + h.c. 



(8) 



The new effective couplings C's contain the information about the underlying theory as 
well as the form factors. The vector couplings are: 



a 



ip) 

VR 



a 



VR 



VL 

(n) 



VL 



q=u,d,s 



q=u,d,s 



q=u,d,s 



q=u,d,s 



il) fil) 

VR JVp 


(9) 


{<}) fig) 

VR JVn 


(10) 


<(<?) fil) 
VL JVp 


(11) 


<(<?) /(<?) 
VL JVn ' 


(12) 
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while the scalar ones read: 



Cgqr (1 - ^ fsp] 

q=u,d,s 

Cgqr (1 - ^ fsl] 

q=u,d,s 

Cgql (1 - fsp) 

q=u,d,s 

Cgql (1 - ^ Z^^^) 

q=u,d,s q=u,d,s 



'-'SR 


= E 

q=u,d,s 




f(g) 


/=<(") 


= E 

q=u,d,s 




Aq) 

JSn 


r'(p) 


= E 

q=u,d,s 




Aq) 

J Sp 




= E 




Aq) 

J Sn 



(13) 
(14) 
(15) 
(16) 



2.2 Transition rates 



The nucleon-level effective Lagrangian can be used to take matrix elements at the atomic 
and nuclear level. In the non-relativistic approximation, the relevant matrix elements are 

{A,Z\^^tlj^\A,Z) = {A-Z)p^''^ 

{A,Z\i;p^%\A,Z) = Zp^P^ 

{A,Z\^P^^'^P^\A,Z) = (A-Z)p(") 

{A,Z\^PN7'^N\A,Z) = 0. (17) 

Here \A, Z) represents the nuclear ground state, with A and Z the mass and atomic number 
of the isotope, while p^^ and p*^") are the proton and neutron densities respectively. The 
conversion rate of the process is written as 

Tconv = ^\CdrD + AGprn^ {rnpC'ilS^P^ + m„4"^ S^^^) + C^^l W^^^ + C^^l 4\/(") " 
+ ^ 



4A4 



Cdl D + AGprn, (m.cf^S^^^ + rrinCfl ^(")) + C^^l AV^^^ + C^^ 4\/(") 



in terms of the dimensionless integrals V'^^\ S^^\ representing the overlap of electron 
and muon wavefunctions weighted by appropriate combinations of protons and neutron 
densities [B]. For phenomenological applications, it is useful to normalize the conversion 
rate to the muon capture rate, introducing the quantity: 

Finally, we note here the branching ratio for the purely radiative process /i — *■ 67 in 
terms of the effective couplings defined above: 

_ r(/i^e7) _ 487r^ , ^ 2. 
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2.3 Sources of uncertainty 

There are two sources of uncertainty in the calculation of the transition rate: (i) scalar 
form factors and (ii) neutron density (for high Z nuclei). The latter uncertainty has been 
carefully discussed in Ref. [6], where several approaches to determine the neutron density 
have been reviewed and used in the calculation of the overlap integrals. Whenever data 
from polarized proton scattering exists, the uncertainty on the overlap integrals 5*'-"^ and 
can be reduced to a few percent even for heavy nuclei such as Pb. Otherwise, it should 
be considered to be of the order of 10%. In this work we focus on the uncertainty induced 
by the scalar density matrix elements in the nucleon. 

The scalar form factors defined in Eq. ([5]) can be re-expressed in terms of ratio of quark 
masses and ratios of nucleon matrix elements as follows ITT]: 



= ^^(1 + 0^ (21) 



ffi = -^{1-0^ (22) 

fsi = ^^y^ (23) 
^ mu + rrid nip 

f'sl = ^^(1-0^ (24) 
niu + rrid nip 

f'sH = -^{l+O'-^ (25) 

niu + nid nip 

fitl = -^y'-^, (26) 

mu + nid nip 



where 



o-nN = ^ {p\uu + dd\p) (27) 

^ ^ {p\uu- dd\p) ^28) 
{p\uu + dd\p) 

y = _J{p\ss\p^_ ^29) 
{p\uu + dd\p) 

Information on the above matrix elements can be obtained from vrA^ scattering data, 
from an analysis of the octet baryon masses within Heavy Baryon Chiral Perturbation 
Theory, or from Lattice QCD. 

For the cr-term, we will use the lattice result 

a^N = (53 ± 2(stat)lf (syst)) MeV , (30) 

whose uncertainty covers determinations from nN scattering |T3|,[T3], from ChPT analysis 
of baryon masses [15], as well as from previous lattice analyses [1^. For the ratio measuring 
isosp in-breaking, we will use [HI |T7j : 

^ = 0.132 ±0.035 . (31) 
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For the ratio y quantifying the strange quark content of the nucleon, the situation is less 
clear. A Chiral Perturbation Theory analysis gives the range y = 0.21 ±0.2 [15]. The large 
uncertainty reflects the poor knowledge of the relevant low-energy constants, even within 
resonance saturation (the matching renormalization scale is arbitrary). A recent lattice 
QCD analysis [12] of the matrix element (A^|ss|A^) within the overlap fermion formulation 
with two dynamical flavors leads to 

y = 0.030 ± 0.016(stat)l[J:J][J^(extrap)t[!:[!J]^(m,) . (32) 

This result is obtained from the lattice matrix element (A^|ss|A^) by dividing out the sigma- 

(s) 

term as calculated in the same lattice simulation. Therefore, the uncertainty in /^^ oc y x 
Utttv is controlled by Eq. [321 with a^^^ = 53 MeV simply providing the normalization. The 
lattice result is consistent with the Chiral Perturbation Theory range, although suggesting a 
much smaller strange content of the nucleon. The difference with respect to previous lattice 
results has been attributed to a lattice artifact (mixing with wrong chirality operator) in 
the Wilson fermion approach. The uncertainty on this value is at the moment dominated 
by statistics. 

For the purpose of this work, we will vary the parameter y within both a "conservative" 
range and an "optimistic" range. For the conservative range we take y G [0,0.4], which 
coincides with the ChPT range of Ref [15] . For the optimistic range we take y G [0,0.05] 
which reflects more closely the recent JLQCD result [12] and seems a realistic guess of the 
uncertainty that will be reached by lattice calculations in the next decade. 

Finally, for the ratios of quark masses, we use the the input [18] 

0.553 ±0.043 (33) 
18.9 ±0.8. (34) 

3 Testing the single operator dominance hypothesis 

We now turn to illustrate the model discriminating power of a combined phenomenological 
analysis of /x — 67 and ^ ^ e conversion on different target nuclei. In order to organize 
the discussion, we deflne here four classes of models, in which only one underlying short 
distance operator dominates over all the others. We call these four classes of models the 
"single-operator" dominance models. We will flrst analyze this simplest class of models and 
then consider the more involved case in which two operators have comparable strengths 
and interference effects cannot be neglected. 

3.1 Dipole, Vector and Scalar models 
• Dipole model 

The Dipole model is deflned by the assumption that, among all LFV short- distance 
operators, the dipole operator is the dominant one. For simplicity, we focus on the 
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case in which the outgoing lepton has definite chirahty [§. Exphcitly, in terms of the 
effective couphngs defined in Eq. [H this class of models is defined by: 

Cd = CdR^Q Celse = Q ■ (35) 

Most supersymmetric scenarios, including SUSY-GUT models [19] and SUSY see-saw 
models [20] fall in this class of models. 



Vector model 1: V^'^'^^ 

This model is defined by the assumption that the transition charge radius operator 
gives the dominant contribution to the LEV lagrangian. The model is defined by 

Cv = Cy]^ = —2 Cyj^ 7^ 0, Celse = , (36) 

and is explicitly realized in large regions of the Left-Right symmetric model parameter 
space [21j. In this model C^^^ ^ 0, while C^yl = 0. 

Vector model 2: V^^^ 

The Vector model 2 is defined by the assumption that the underlying dominant 
operator is an effective Z-penguin. The ratios of couplings of different quarks is 
governed by the couplings of the Z° coupling to quarks. The model is defined by: 

Cv ^ 4t = ^ ^ 0, C,,.e = , (37) 
a 

where a is the ratio of the down and up quarks coupling to the Z-boson: 

Tl+T^^-{Qu, + Qu,)sm'ew 

With this value of a (corresponding to sin^ 9w = 0.223) we obtain Cy]^/Cy]^ = —9.26, 
in contrast to the V^'^^ model. 

Scalar model 

This model is defined by: 

Cs = C^SR ^ ^SR ^ ^fli 7^ 0' Celse = . (39) 

This model may be explicitly realized in some regions of the usual R-parity conserving 
SUSY see-saw parameter space [22] (large tan/? and relatively low "heavy" Higgs 
sector) and within R-parity violating SUSY 



^Allowing for the presence of outgoing leptons with both chirahties (e.g. both Cdr 7^ and Cdl 7^ 0) 
would not change the conclusions of the single-operator analysis of this section. 
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Figure 1: The y-parameter dependence of the conversion branching ratio in the scalar 
model. 

Among the above models, the scalar model suffers from the uncertainty in the y pa- 
rameter. We show in Fig. [1] the y-parameter dependence of the conversion branching ratio. 
The uncertainty is quite large if we take the conservative range, y G [0, 0.4]. 

Each of the above classes of models has only one free parameter - the ratio Ci/ h? of 
the dominant effective coupling over the square of the new physics scale. It is clear, then, 
that the single-operator dominance hypothesis makes parameter-free predictions for ratios 
of LFV branching fractions and therefore it can be tested so long as two LFV rates are 
measured. We will discuss how well one can distinguish models in the presence of the 
theoretical uncertainties. 



3.2 n ^ e'j n ^ e conversion 

If yU — > 67 and fi —>■ e conversion in at least one target nucleus are observed, this immedi- 
ately opens up the possibility to test the Dipole dominance model. In fact, in this model 
the ratio 

R(^Z) = ^Jj-lA^ (40) 

is entirely fixed by the overlap integrals D |6], which are essentially free of theoretical 
uncertainty. R{Z) is predicted to scale as C(a/7r) and we plot it in Fig. [2l We omit 
from the plot the points corresponding to es^Er, ys^Ta, and 79^ Au, as data on the nucleon 
densities are either obtained from quite old experiments or not well established [2lj. Any 
deviation from the pattern shown in Fig. [2] would imply the presence of scalar or/and 
vector contributions. In order to disentangle these possibilities, one needs to study the 
target dependence of the conversion rate. 
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Figure 2: Ratio R{Z) of /x 
dominance model. 



e conversion over — e'j) versus Z in the case of Dipole 



3.3 Target dependence of /i ^ e conversion 

In principle, any single-operator model can be tested with two conversion rates, even if 
/i — > 67 is not observed. To illustrate this point, we update the analysis of Ref. [6] and 
plot in Fig. [3] the conversion rate (normalized to the rate in Aluminum) as a function of 
the Z of the target nucleus, for the four classes of single-operator models defined above. 
Compared to Ref. [6], the novelty here is the inclusion of a second vector model {V^^'^). 

The results of Fig. [3] show some noteworthy features. First, we note the quite different 
target dependence of the conversion rate in the two vector models considered. This can be 
understood as follows: in the case of the V^"'^ model, the behavior in Fig. [3] simply traces 
the Z-dependence of V^^^ (the photon only couples to the protons in the nucleus). On 
the other hand, in the case of the V^^^ model, the Z boson couples predominantly to the 
neutrons in the nucleus and the target dependence of the ratio V^'^^V^^^ ~ (A — Z)/Z 
generates the behavior observed in Fig. [31 

Next, let us focus on the actual discriminating power of the Z-dependence. Clearly, 
the plot shows that the model-discriminating power tends to increase with Z. This is 
a simple reflection of the fact that the whole effect is of relativistic origin and increases 
in heavy nuclei. So in an ideal world, in order to maximize the chance to discriminate 
among underlying models, one would like to measure the conversion rate in a light nucleus, 
say Aluminum or Titanium, as well as in a large-Z nucleus, like Lead or Gold. This 
simplified view, however, has to be confronted both with theoretical uncertainties and the 
actual experimental feasibility. Concerning the uncertainties, a simple analysis shows that 
the dominant uncertainty coming from the scalar matrix elements almost entirely cancels 
when taking ratios of conversion rates (even using the conservative range y G [0, 0.4] for 
the strange scalar density matrix element). Moreover, in the large-Z tail of the plot, some 
residual uncertainty arises from the input on the neutron density profile. When polarized 
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z 

Figure 3: Target dependence of the /i — >■ e conversion rate in different single-operator 
dominance models. We plot the conversion rates normalized to the rate in Aluminum 
{Z = 13) versus the atomic number Z for the four theoretical models described in the 
text: D (blue), S (red), V^^^ (magenta), V^^^ ( green). The vertical lines correspond to 
Z = 13 (Al), Z = 22 (Ti), and Z = 83 (Pb). 

proton scattering data exists, the uncertainty on the ratios of conversion rates becomes 
negligible. This point is illustrated by Tabled], where we report the detailed breakdown of 
uncertainties in the ratios -B^^e(Ti)/-B^-^e(Al) and i?^_>e(Pb)/-B^^e(Al). For other targets, 
the uncertainty induced by neutron densities never exceeds 5% i6j. The conclusions of this 
exercise are that: 

• The theoretical uncertainties (scalar matrix elements and neutron densities) largely 
cancel when we take a ratio. 

• As evident from Fig. [3l a realistic discrimination among models requires a measure 
of i?^^e(Ti)/i?^_^e(Al) at the level of 5% or better, or alternatively a measure of 
-B^^e(Pb)/i?^_>e(Al) at the 20% level. These are two cases that well represent the 
trend in light and heavy target nuclei. 
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S 


D 






B(At^e,Ti) 


1.70 ± O.OOSy 


1.55 


1.65 


2.0 


B{p^e,Al) 


0.69 ± 0.02p„ 


1.04 


1.41 


2.67 ± 0.06p„ 



Table 1: Ratios of conversion rates in Titanium and Lead over Aluminum, in each of 
the four single-operator models: scalar (S), dipole (D), vector 1 (photon coupling to the 
quarks) and vector 2 {Z boson coupling to the quarks). In the scalar model, the scalar 
form factor induces a negligible uncertainty in the ratios involving two targets (denoted by 
the subscript y) . In the case of Lead over Aluminum, the small uncertainty is dominated 
by the neutron density input (denoted by the subscript p„). 

4 Testing the two-operator dominance hypothesis 

In the last section we have discussed how to test the hypothesis of a single operator dom- 
inance, and how to discriminate among different single-operator dominance models. If 
the single operator dominance hypothesis fails, one is lead to consider next simplest case, 
namely the two-operator dominance models, defined by the assumption that only two 
underlying operators have appreciable coefficients. Each model is characterized by two 
parameters, the effective strength Ci/A^ of one of the two operators and the ratio C2/C1 
of the effective couplings of the two dominant operators. This class of models can be tested 
so long as two double ratios of LFV rates are available (three LFV measurements!). 

For the sake of illustration, we will consider the following three two-operator models: 
Dipole-Scalar, Dipole- Vector (Z) and Scalar- Vector (Z). We consider both the case of con- 
structive and destructive interference among the two dominant operators, assuming that 
the ratio of Wilson coefficients r = C2/C1 is real (a relative phase can be included but 
it would unnecessarily complicate the analysis at this early stage). In order to test this 
class of models, one has to assume that at least three LFV processes have been observed, 
so one can construct two independent double ratios that are entirely determined by the 
single parameter r. In models involving the Dipole operator among the dominant terms 
(such as Dipole-Scalar and Dipole- Vector) the three observables could be (i) yU ^ 67 and 
fi —>■ e conversion in two different targets; or (ii) fi —>■ e conversion in three different targets. 
In models that do not involve a Dipole term (such as Scalar- Vector), only the possibility 
(ii) above is available. As representative target nuclei, we have chosen aluminum (Al), 
titanium (Ti), and lead (Pb). 

4.1 Dipole-Scalar 

In terms of the parameters introduced in Section 13.11 this model is defined by C5 7^ 
and Cf) = ±^ Cs- The single-operator models are recovered in the limiting cases r — 
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(scalar) and r oo (dipole) i Note that in this particular case the asymptotic dipole 
regime is reached already for r ^ 1 because of the peculiar normalization of the scalar 
operators (suppressed by the factor Gpfngm^). 




Log r Log r 



(b) (c) 

Figure 4: Dipole-Scalar model: Ratios -B^^e(Al)/i?^^e7 (panel (a)), i?^^e(Pb)/i?^^e(Al) 
(panel (b)), and -B^^e(Ti)/-B^_+e(Al) (panel (c)) as a function of Log]^Q(r) for positive 
Cd/Cs- See text for details. 

We illustrate the features of this model in Figs. H] and [5l which correspond to pos- 
itive and negative sign of the ratio Cn/Cs, respectively. Panel (a) shows the behav- 
ior of i?^^e(Al)/B^^e7 versus the parameter r, while panels (b) and (c) show the ratios 
-B^^e(Pb)/i?^^e(Al) and -B^^e(Ti)/_B^^e(Al), respectively. In panels (a) and (c) the curve 
is widened in the interference region by the uncertainty in the scalar form factors. The 
dominant uncertainty comes from the input parameter y, characterizing the strangeness 
content of the nucleon. On the other hand, the ratio i?^^e(Pb)/i?^^e(Al) is affected not 

•^We consider here the case in which dipole and scalar operators produce outgoing lepton with definite 
chirality (L or R). If both chiralities are allowed, then in principle Cdr/Csr ^ Cdl/Csl and one more 
parameter has to be introduced in the analysis. 
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only by the uncertainty in the scalar form factors, but also by the uncertainty induced in 
the overlap integral by the neutron density in Pb. The width of the bands in panel (b) is 
determined by the most conservative combination of two kind of uncertainties. 

In all panels the wide band corresponds to the range y G [0,0.4], while the narrow 
band corresponds to the range y G [0,0.05]. This illustrates the effect of current and 
future hadronic uncertainties on the process of extracting information on short distance 
LFV couplings. The prominent feature in Fig. [5] is induced by the destructive interference 
dipole and scalar amplitudes. 




(b) (c) 

Figure 5: Dipole- Scalar model: Ratios i?^^e(Al)/-B^^e7 (panel (a)), -B^^e(Pb)/i?;i^e(Al) 
(panel (b)), and i?^^e(Ti)/i?^^e(Al) (panel (c)) as a function of Logio(?") ^r negative 
Cd/Cs- See text for details. 

4.2 Dipole- Vector 

In terms of the parameters defined in Section 13. this model is defined by Cy 7^ and 
C/) = Cy- The single-operator models are recovered in the limiting cases r ^ 
(vector) and r — > 00 (dipole). In figures E] and [7] we plot the ratios -B^^e(Al)/i?^^e7 
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(panel (a)), S^^e(Pb)/5;,-.e(Al) (panel (b)), and B^^e(Ti)/5^^e(Al) (panel (c)) versus 
the parameter r. Figures [6] and [7] correspond to positive and negative sign of the ratio 
Cd/Cvi respectively. Within this model, the only source of uncertainty arises from the 
vector overlap integral l^*-""-' (Pb) , sensitive to the neutron density in Pb. This uncertainty 
is quantified by the thickness of the band in panel (b). 




-2 -1 1 2 3 4 -2 -1 1 2 3 4 



Log r Log r 

(b) (c) 

Figure 6: Dipole- Vector model: Ratios i?^^e(Al)/-B^^e7 (panel (a)), i?^^e(Pb)/5^^e(Al) 
(panel (b)), and i?^^e(Ti)/i?^_»e(Al) (panel (c)) as a function of Logxo('") positive 
Cd/Cv- See text for details. 
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(a) 




Logr 



Logi- 



(b) (c) 

Figure 7: Dipole- Vector model: Ratios -B^_,e(Al)/i?^_»e^ (panel (a)), i?^^e(Pb)/i?^^e(Al) 
(panel (b)), and i?^^e(Ti)/i?^_+e(Al) (panel (c)) as a function of Logio(^) negative 
Cd/Cv- See text for details. 



4.3 Scalar- Vector 

In terms of the parameters defined in Section I3.H this model is defined by Cy 7^ and 
Cs = ±r Cy- The single-operator models are recovered in the limiting cases r (vector) 
and r — > 00 (scalar). Since the Dipole term is assumed to be subdominant, in this case 
we include in the analysis only the ratios -B^^e(Pb)/_B^^e(Al) and i?^^e(Ti)/-B^_^e(Al), 
shown in panels (b) and (c) of Figures [8] and [9] (for positive and negative values of Cs/Cy, 
respectively). While the ratio B^^e(Ti)/i?^^e(Al) is affected only by the uncertainty in y, 
the ratio -B^^e(Pb)/-B^^e(Al) is affected also by the uncertainty in the Pb neutron density 
(through the overlap integrals). The width of the bands in the plots is determined by the 
most conservative combination of two kind of uncertainties. 

In all panels the wide band corresponds to the range y G [0,0.4], while the narrow 
band corresponds to the range y e [0, 0.05]. As in the case of the Dipole-Scalar model, the 
bands illustrate the effect of current and future hadronic uncertainties on extracting short 
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distance LFV couplings. 




Log r 



Logr 



(b) 



(c) 



Figure 8: Scalar-Vector model: Ratios -B^^e(Pb)/-B;x^e(Al) (panel (b)) and 
S^^e(Ti)/i?^^e(Al) (panel (c)) as a function of Log]^Q(r) for positive Cs/Cy- See text 
for details. 



T 



7 

Log r 




(b) 



(c) 



Figure 9: 



Scalar-Vector model: Ratios i?^^e(Pb)/i?^^e(Al) (panel (b)) and 



i?^^e(Ti)/i?^^e(Al) (panel (c)) as a function of Logio(^) negative Cs/Cy. See text 
for details. 



We conclude this section by summarizing what one could learn about the two-operator 
dominance models in the case that two double ratios of LFV rates could be measured 
experimentally. Our exercise shows that: 

• The current theoretical uncertainty on the strange content of the nucleon prevents a 
realistic test of the two-operator models involving the Scalar amplitude. The range 
y G [0, 0.4] induces uncertainties of up to one order of magnitude in the relevant 
double ratios in the interference region (thick bands in all plots above). However, 
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Figure 10: The pseudoscalar-Higgs mass dependence of i?^^e(Al) / -B^^e7 in a SUSY model 
with a left-handed slepton mixing. The slepton masses, the gaugino masses at the GUT 
scale, and the Higgsino mass parameter are all fixed to be 1 TeV. The light and dark 
shaded regions respectively correspond to the conservative and optimistic ranges of the y 
parameter. 



the uncertainty within reach of lattice QCD calculations will remove this obstacle in 
the coming years (this is illustrated by the thin bands in all plots above). 

• Testing and discriminating among two-operator dominance models requires an ex- 
perimental precision on the LFV rates that is comparable to the one needed to test 
the single operator models. 

5 Application to a SUSY model 

An example of the two-operator dominance model is given by a SUSY scenario with flavor 
mixing in the left-handed sleptons. Such a mixing, for example, can be induced from the 
Yukawa interaction in the see-saw model. As it is shown in Ref. [22], the scalar operator 
originated from the Higgs-boson-exchange diagrams can be sizable in this model if tan/3 
is large and the heavy Higgs boson is relatively lighter than the other SUSY particles. 
The ratio i?^^e(Al)/-B^_te7 can therefore be enhanced in such a parameter region, while 
the ratio i?^^e(Pb)/-B^^e(Al) can show substantial deviations from the dip ole- dominance 
value. 

In Fig. [To] we show the pseudoscalar-Higgs mass (m^) dependence of the ratio for 
/i > (left) and < (right). We have taken the common mass (msusY = 1 TeV) for 
the slepton masses, the universal gaugino mass at the GUT scale, and the Higgsino mass 
parameter, and we fixed tan/3 = 60. Since the scalar operator does not decouple in the 
'^susY oo limit, we see the enhancement in the small tua region. The light (dark) shaded 
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Figure 11: The pseudoscalar-Higgs mass dependence of B^^e{Ph)/ B^^ei-^^) in a SUSY 
model with a left-handed slepton mixing. The slepton masses, the gaugino masses at the 
GUT scale, and the Higgsino mass parameter are all fixed to be 1 TeV. The light and dark 
shaded regions respectively correspond to the conservative and optimistic ranges of the y 
parameter. 



regions correspond to the conservative (optimistic) range of the y parameter, y G [0, 0.4] 
{y G [0,0.05]). Within the same framework, the ratio B^^e{Ph)/B^^f;{Al) is shown in 

In both cases, the theoretical uncertainty becomes significant as the scalar operator gets 
important. In the context of this explicit supersymmetric model, a precise determination 
of the y parameter is quite important in order to extract information on the underlying 
model parameters. To illustrate this even more explicitly, in Fig. [T2] we show for fi > 
the niA dependence of i?^^e(Al)/5^^e7 (left) and i?^^e(Pb)/-B;i^e(Al) (right) for different 
values of tan/3 = 40, 50, 60. In these plots, only the small uncertainty window is reported 
{y G [0,0.05]), to illustrate the enhanced discriminating power. 



6 Conclusions 

In this work we have investigated whether the target- dependence of /i-to-e conversion 
rate can be exploited to discriminate among underlying dynamical mechanisms of lepton 
flavor violation, once one takes into account realistic hadronic and nuclear uncertainties. 
The major source of theoretical uncertainty arises from the nucleon matrix element of the 
strange quark scalar density. This is expressed in terms of the parameter y (see Eq. I29l) . 
which we have varied within two ranges reflecting the current uncertainty {y G [0, 0.4]) and 
the projected uncertainty within reach of lattice QCD calculations {y G [0,0.05]). 

In order to assess the model discriminating power of a combined phenomenological 
analysis of —>■ and —>■ e conversion on different target nuclei, we have defined four 
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Figure 12: The pseudoscalar-Higgs mass dependence of i?^_^e(Al)/_B^^e7 (left) and 
-BM^e(Pb)/-B^^e(Al) (right) for different values of tan/3. The thickness of the bands cor- 
responding to varying y in the optimistic range [0, 0.05]. 



classes of models, in which only one underlying short distance operator dominates over all 
the others (Dipole, Scalar, Vector (7) and Vector {Z)). Ratios of LFV branching fractions 
can be used to test the various models. The single-operator hypothesis can be tested with 
at least one ratio (two LFV measurements), while the two-operator models, where two 
operators have comparable strength and interfere, can be tested with at least two ratios 
(three LFV measurements). 

Our conclusions are encouraging: the theoretical uncertainties (even at the current 
level) are not an issue in testing the single-operator dominance model, as they largely 
cancel when we take ratios of different conversion rates. On the other hand, the current 
uncertainty prevents meaningful tests of two-operator models involving the Scalar operator, 
as it produces errors of up to one order of magnitude in the double ratios in the interference 
region. However, with the anticipated reduced lattice error on the strange content of the 
nucleon, this will not be an issue in the future. We have illustrated these main conclusions 
also in the context of a suprsymmetric model. 

Having established that the hadronic uncertainties will not be a limiting factor, we can 
ask how well one should measure the LFV rates in order to discriminate the underlying 
models. Fig. [3] shows that a realistic discrimination among single-operator models requires 
a measure of the ratio of conversion rates in light nuclei (such as i?^^e(Ti)/-B^^e(Al)) at 
the level of 5% or better. Alternatively, one would need to measure the ratio of conversion 
rates in a heavy and light element (such as i?^^e(Pb)/i?^^e(Al)) at the 20% level. Similar 
accuracy is required on the experimental side to be sensitive to interference effects when 
more than one operator is at work. Whether these challenging benchmark numbers can be 
reached in the future round of experiments [251 126] depends on many issues, including the 
value of the branching fraction themselves (we are concerned here with ratios). Nonetheless, 
we hope that our results will stimulate further experimental efforts towards measurements 
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of /i ^ 67 and /x-to-e conversion and consideration of various options for target nuclei. 
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